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We study dc conductivity of a Weyl semimetal with uniaxial anisotropy (Fermi velocity ratio 
^ = v±/v\\ 7^ 1) considering the scattering of charge carriers by a wide class of impurity potentials, 
both short- and long-range. We obtain the ratio of transverse and longitudinal (with respect to the 
anisotropy axis) conductivities as a function of both ^ and temperature. We find that the transverse 
and longitudinal conductivities exhibit different temperature dependence in the case of short-range 
disorder. For general long-range disorder, the temperature dependence (~ T^) of the conductivity 
turns out to be insensitive of the anisotropy in the limits of strong (^ ^ and ^1) and weak (^ ?^ 1) 
anisotropy. 
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I. INTRODUCTION 

Weyl semimetals (WSM) are three-dimensional (3D) 
analogs of grapheneii^ Their quasiparticles are described 
by the massless 3D Dirac Hamiltonian. Such systems 
were first proposed as an exotic theoretical possibil¬ 
ity, with an expectation to observe some of their fea¬ 
tures in pyrochlore iridates^ and in certain semiconduc¬ 
tor heterostructuresRecent experiments uncovered 
several chemical compounds, which can be classified as 
WSM^— . The subject is attracting growing theoreti¬ 
cal and experimental interest, (see, e.g., the reviews in 
Refs, [lllilii). 

Charge transport in WSM has also received consid¬ 
erable attention and a number of interesting phenomena 
have been discovered (see Ref. [13 and references therein). 
The main emphasis was on manifestations of topological 
effects: an additional topological protection of the gap¬ 
less spectrum near the Dirac points in 3D and the chi¬ 
ral anomaly^”— as well as an unusual Kondo effecfi^^. 
Disorder and impurity effects have also attracted signif¬ 
icant attentionj22ri2& A detailed study of the influence of 
Coulomb disorder in highly compensated WSM was un¬ 
dertaken in Ref. [l^. The problem is that even in the 
simplest case of delta-correlated disorder, the field theory 
of a WSM is non-renormalizable. Progress in overcor ning 
this stumbling block was recently achieved in Refs. 1^1^ . 
where the conductivity of weakly-disordered semimet¬ 
als was treated in terms of the ^-expansion within the 
renormalization group (RG) approach. However, Weyl 
semimetals correspond to 5 = —1 and the predictions 
are qualitative. Nevertheless, it appeared to be possible 
to reveal some specific features of the conductivity for 
different kinds of disorder and even to predict a disorder- 
driven quantum phase transition distinct from the Ander¬ 
son transition. A RG approaches has been also applied to 
demonstrate the possibility of non-Fermi-liquid behavior 
in systems with isotropic 3D Dirac points having a weak 
disorder and Coulomb-type interparticle interactions 


The computation of the isotropic conductivity was re¬ 
cently performedeSiSiieS in the framework of the kinetic- 
equation technique, improved by the self-consistent Born 
approximation, as well as numerically^SS 

Most Weyl semimetals obtained in laboratories so far 
are anisotropicThe anisotropy can be also induced 
by a linearly-polarized electromagnetic wave^^ Some ef¬ 
fects related to the anisotropy and tilting of the Dirac 
cones have recently been treatedHowever, the in¬ 
fluence of anisotropy on the charge transport in WSM 
with quenched disorder has remained unexplored so far. 
In this paper, we analyze the effects of disorder and 
anisotropy on the conductivity, in the framework of the 
Born approximation. We treat these effects in the dia¬ 
grammatic approach, both for short- and long-range dis¬ 
order taking into account the uniaxial anisotropy. We 
obtain analytical results for the cases of strong and weak 
anisotropy. The diagrammatic framework allows us to 
find out the validity range of the Fermi-liquid approach 
in WSM. 

The potential produced by impurities (disorder poten¬ 
tial) is introduced as a rather general spatial distribu¬ 
tion u{t). Its characteristic scale in momentum space 
is specified by a parameter po. Thus, at small values 
of Po, we are dealing with long-range disorder, whereas 
large Po corresponds to short-range disorder. The ratio of 
Fermi velocities perpendicular and parallel to the z axis, 
^ = v±/v\\, becomes an additional control parameter of 
the problem. 

We have found that for sufficiently low temperatures 
(or doping level 5 f), when the disorder can be consid¬ 
ered to be a short-range one (po ^ max{T/i;||,5F/'^}), 
the longitudinal and transverse conductivities exhibit a 
different temperature dependence 

15 ('^llPo)^ 

where p(p^/Po) is the impurity structure factor and 
dg{0) = dg{x)/dx\x=o- The temperature dependence in 
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Eq. (dJ) is in fact the first term of the asymptotic series 
in max{T/i;||,5F/'^}/po ^ 1- This dependence becomes 
even more pronounced at higher temperatures T ^ pqV\\ . 

The ratio saturates to a temperature- 

independent constant in the high-temperature (or long- 
range disorder) limit, when (po ^ max{T/'r’||,At 
5 f ^ T, the temperature dependence of both compo¬ 
nents of the conductivity obeys the relation 

(T||, cr_L ~ (2) 


potentials. We compute the conductivity for different 
T/sy ratios and different values of the anisotropy pa¬ 
rameter The conclusions are presented in Section V. 
Several important technical issues are presented in the 
Appendices. 

II. FORMALISM 
A. The model 


regardless of the particular form of the disorder potential 
or impurity structure factor. In the long-range disorder 
limit, for ^ ^ 1 , ~ 1 , and ^ ^ 1 , we found that the 

ratio approaches constant values, independent 
of temperature and disorder potential. The results con¬ 
cerning the conductivity ratio can be summarized in the 
form 


^.2 r3/(i + 41n2)+0(r'), e»l, 

^ 1 + e = i + <5^, (3) 

Ic 3 + (!?(e 2 ), ^< 1 , 


where C 3 is a constant of the order of unity. 

The scales of the problem are the temperature T of the 
system, Fermi energy 5 f, and the elastic scattering rate 
l/r. The latter is assumed to be small, so 


maxj^F, T}r 

is a small parameter of the problem. Our task is to com¬ 
pute the Drude conductivity. The fermion doubling the¬ 
orem implies that the WSM spectrum always has an even 
number of Weyl points. For simplicity, we assume that 
the charge carriers have identical spectra near all Weyl 
points and compute the conductivity per point All the 
results for the conductivity should be multiplied by the 
number of Weyl points. 

The disorder potential is assumed to be quite weak, so 
that the Born approximation is justified 


|w(r)| < inin{?;x,V||}po. (5) 

Note here that our results are applicable not only to the 
Weyl semimetals themselves, but to a wider class of ma¬ 
terials with the 3D Dirac spectrum (topological Dirac 
semimetals etc., see the classification given in Ref. [3^ . 

In Section II, we describe the formalism used in dealing 
with the scattering problem, give all the necessary defini¬ 
tions, and describe the employed calculation procedures. 
In this section, we also analyze the limits of applicability 
of the Fermi-liquid approach in the case of WSM. In Sec¬ 
tion III, we revise the general temperature dependence 
of the conductivity in isotropic Weyl semimetals with 
long- and short-range disorder. Section IV deals with 
the anisotropic case. We provide a detailed analysis of 
the solution to the Dyson equation for the singular part 
of the vertex function for short- and long-range disorder 


We consider a system with a Hamiltonian, which gen¬ 
erally has uniaxial anisotropy. Namely, the Fermi ve¬ 
locity along a specified axis no (the longitudinal compo¬ 
nent) is different from that in the perpendicular direction 
v± = Ci;|| 


H — Hq + i^dis; 

Ho = -i j V’^(r)[v||(T||5r|| +v:i_(T^^r^]'lp{v)dY, 

Hdis = J tp\r)u{r)'tp{r)dr, 

where u{t) is the disorder potential, are quasiparti¬ 
cle field operators, and = crno, ry = rno, cr± = 
a — trno, = r — no (nor) are the projections of Pauli 
matrices and radius vectors. The disorder potential cor¬ 
relation function reads 

[ dre-^^’^{u{r)u{0)) = ^^5 (g) , (7) 

J Po \Po/ 

where g{p‘^IP q) is the dimensionless Fourier transform 
of the normalized (^(0) = 1) disorder structure factor. 
In our case, it incorporates the correlation of impurity 
positions, as well as the form of the potential. The pa¬ 
rameter 1^0 plays the role of the amplitude of the disor¬ 
der potential and nimp is the concentration of impurities. 
The potential is assumed to be isotropic. The last as¬ 
sumption certainly depends on the nature of the disor¬ 
der. For example, the isotropic Coulomb potential may 
acquire an anisotropic screening by charged carriers with 
an anisotropic spectrum. Here, we neglect the induced 
anisotropy of the potential. Our task is to emphasize the 
main effect of the anisotropic spectrum on the transport 
properties, which manifests itself even with an isotropic 
potential. 

To simplify the analysis, we rescale the coordinates and 
pj operators to absorb the anisotropy in Hq according to 


r-ii = rfi, rj_ = ip{r) = 

Then, the Hamiltonian changes its form to 

Hq = -W|| j ij'\r'){cTd'^)i;'{r')dr', 

H'dis = J drV't(r')M(r|| + {r '), 


(8) 


( 9 ) 
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which leads to the modified disorder correlation function 


J dr'e 'P--' (w(r' +^rY)w(0)) = g' 
fp^\ _ /(i-r")(pno)" + rv\ 

^Wo) H J’ 


( 10 ) 


where we have decomposed the momentum as: ~ 

(p2-pp->(pno)^+r^[p^-(pno)^]. 

The Feynman rules are extracted from m and (iini, 
and depicted in Fig. [T] 


e,q 



1 


£ - n||pcr + iO 

^imp^o f (1 -C~^)(pno)^ 
Po 


(a) G(£,q) = —= 

q-p 

/ \ 

(b) E(£,q)= ^ 

Figure 3: (a) Dyson equation for the exact fermion 
Green’s function; (b) Fermion self-energy in the 
one-loop approximation. 




Figure 4: Dyson equation for the vertex 
Ja{e,e 


Figure 1: Feynman diagrams for the non-interacting 
fermion retarded Green’s function and the disorder 
correlation function. 


B. Conductivity 


The conductivity tensor is found via the Kubo formula 

( 11 ) 


O’a/3(w,0) = 


lUJ 


n 


(w) =i jdtdr{[ja{t,r),ji 3 {O,O)])e^^* 0 {t), 


where (cj) is the retarded polarization operator and 
ja{t,r) = V’l(l,r) 1*) is the quasiparticle current 

operator. The averaging (...) is assumed to be done over 
the Gibbs distribution as well as over different realiza¬ 
tions of disorder. We will be interested in the system re¬ 
sponse to a uniform electric field, constant in time. For 
this purpose, we set q = 0 and cj -> 0. 






Ga{e,ci) = 


exact Green’s function 


The computation of the polarization operator in the 
lowest order of the 1 / (sft) expansion involves the sum¬ 
mation of disorder ladder series (accounting for the dif¬ 
ference between averages {jaj^) and {ja){ji3)) as well as 
one-loop corrections to the Green’s functions (mainly re¬ 
sponsible for the finite quasiparticle lifetime). 

Due to the Onsager relations, the conductivity (as 
well as the polarization operator) is a symmetric tensor. 
Thus, the exact diagrammatic representation for 
(we denoted 11 ^ 130 ) = allows for its symmet¬ 

ric form (see Fig. [2]). 

The non-interacting and disorder-averaged Green’s 
functions are related via the standard diagrammatic 
equation depicted in Fig. [3^ and have the form 

G^{s, q) = [e - - S^(£, q)] (12) 

The self-energy is worked out in the Born approximation; 
its diagrammatic representation is shown in Fig. [3]3. 

The ladder summation in the polarization operator is 
included into the renormalized vertex J corresponding to 
the electric current. To perform the summation, one has 
to solve the corresponding Dyson equation depicted in 
Fig. m 

The conductivity tensor in a system with the uniax¬ 
ial anisotropy is characterized by only two eignenvalues 
(T \\, ( 7 ±_ and in “principal rescaled axes” takes the form 

/^V|| 0 o\ 

= 0 (T_l 0 . (13) 

V 0 0 


Jq,(£, £ + u;, u;; p) = 



exact vertex 


The general expression for the conductivity in the leading 
1 / (^F^) approximation reads 


£,p 


Figure 2: Diagrammatic representation of the 
polarization operator. 



ds Q-gff (e) 

27r cosh^ ^ ’ 


where 


( 14 ) 
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(^) 


2 2 
e^vfi 


Att 


dp 

(27r)^ 


tr|G^(e,p) [J^^^(e,p) 


(T„]G^(e,p)(T/3 + G'^(e) [J^^^(e,p) - ap]G^{e,p)aa 
-[G^(e,p) - G^{e,p)]aa [G^(e,p) - G^(e, p)]cr/3|. 


(15) 


Here, for brevity, we denoted 5,0; p) = 

J^^^(5,p), and is the singular part of the vertex 

function (see Appendix A for details). 


C. Validity of the Fermi-liquid approach 

Expression m for (7aj3{^) is rather complicated and 
deserves special attention. The integrand defining (Ta^{s) 
is deliberately split into two lines. The second line leads 
to a convergent integral because (G^ — ^ 1/p^, 

when p ^ oc. The first line of ([T5|) though, contains 
the vertex which makes it problematic even in the 

isotropic case {v\\ = v± = v). At first glance, the expres¬ 
sion seems suitable enough, because the product is 

sharply peaked at p = = s/v. The integral therefore, 

appears to be completely determined by the vicinity of 
p = Ps. However, a more attentive look reveals that the 
product G^G^ converges rather slowly: G^G^ ^ 1/p^ 
at Pe ^ p ^ T^Pe, and G^G^ ^ 1/p^ at p rep^. 
Thus, we need to know at least the asymptotic behavior 
of J^^^{e^p) at large momentum in order to complete 
the calculation. 

This difficulty is closely related to the justification of 
the Fermi-liquid approach. Indeed, this justification in¬ 
cludes two constraints: (a) the smallness of the quasipar¬ 
ticle scattering rate guarantees the applicability of per¬ 
turbation theory (in most cases, the vacuum state, the 
Fermi sphere, survives); and (b) all observables are de¬ 
termined by the scattering processes in the proximity of 
the Fermi sphere. The latter statement allows expanding 
the momentum near the Fermi surface in every integra¬ 
tion (d^p = {e^/v^dp in 3D) making any momentum 
integral one-dimensional and fast converging. 

The justification of the second constraint follows from 
the fact that a momentum integral defining the ob¬ 
servable always involves the G^{p)G^{p) term (sharply 
peaked near the Fermi sphere). However, we have just 
seen that the latter argument fails when computing the 
WSM conductivity. The integrand defining the conduc¬ 
tivity does contain a sharp G^(p)G^(p) peak as well as 
a long tail (see Fig. 

To justify the Fermi-liquid framework, we need to 
accurately estimate the contribution given by the tail. 
We divide the integration domain according to p^dp = 


(p^ — Pe^)dp-h Pe^dp. Then, is split into two parts 



I = y(P -Pe‘^)G^i£,p)[Ja^^{s,p) - aa]G^{e,p)a/3dp, 

II = ps^ j G^{s, p) - aa]G^{e, p)(7p dp. 

(16) 

Part H is already convergent. The convergence of part I 
entirely depends on the asymptotic behavior of the ver¬ 
tex function. The ultraviolet behavior of the vertex is 
discussed in Appendix C, where it is proved to decay fast 
enough to secure its convergence. Part I is estimated in 
Appendix D, where it is shown that 

I-max|l,^| —II. (17) 

I Ps) er 

This means that the Fermi-liquid approach holds with 
l/(max{5F, T}r) accuracy. Therefore, some care needs 
to be taken when working out the weak-localization cor¬ 
rections (WLC) to conductivity in, e.g., the 3D case. 
Indeed, this correction is determined by the next term 
in the 1/(5 ft) expansion, and the WLC in 3D reads 
^o^^^Gc/o^ ^ 1/(5ft). Therefore, when computing cor¬ 
rections to the conductivity in the 3D case, one has to 
take into account the values of the momentum far from 
the Fermi surface. 

As a result, expanding all the momentum integrals near 
the Fermi surface, we can omit (in the leading order) the 



Figure 5: Contributions to the conductivity. Region I 
represents the Fermi liquid contribution; region H 
originates from the scattering far from the Fermi 
surface. 
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and G^G^ terms in the expression (113 for the 
conductivity and obtain the standard formula 

2 2 

+ G^(£)J_g^^(e,p)G^(e,p)c7„|. 

( 18 ) 

III. CONDUCTIVITY IN THE ISOTROPIC 
CASE 

In order to understand how to tackle the anisotropic 
problem, we briefly outline the diagrammatic deriva¬ 
tion of the isotropic conductivity. Wherever isotropy 
is implied, we omit the || and ± indexes, and write 
= v± = V. In the corresponding limits, we reproduce 
the earlier results for the condnctivity j^^i^^i^^ 


A. Disorder averaging 

The self-energy of the Green’s function has matrix 
structure (due to the presence of a-matrices) and consists 
of two parts: -1- The first part is respon¬ 

sible for the renormalization of Fermi velocities and the 
rescaling of the Fermi field operators; the second part is 
responsible for the non-vanishing relaxation rate. Pertur¬ 
bation theory can be applied only if \T^r\ maxj^F, T}. 
The estimate gives 

|sfl fvpo\ 1 w ^ • /i 

—— - j — ^ A(po) inin i 1 

1 (^impMo) min{e2, (pov)^} 

- -—- 

r v^Pq 

where |S^| is understood as the modulus of any of its 
matrix component, 5 ^ max{T/i;, (^f can be con¬ 

sidered here as the the doping level). The dimensionless 
coupling constant 



A(.o) = ^ (20) 

represents the strength of the potential. It is defined at 
the characteristic disorder scale pq. 

Clearly, despite the smallness of the Fermi-liquid pa¬ 
rameter 1 / {sr) , the first part of the self energy may 
lead to a strong renormalization for short-range {po ^ 
s/v) disorder, depending on the initial coupling strength 
A(po)- This is in contrast to the situation in ordinary 
metals, where the ultraviolet cut-off ypo/s-p ^ 1 and the 
smallness of the 1 / (sft) parameter guarantees the valid¬ 
ity of perturbation theory. In WSM, it is appropriate to 
consider both limits: short- (po ^ max{pF, T/i;}) and 
long-range (po ^ max{pF, T/'i;}) disorder. 


1. Renormalization due to short-range disorder 


The analysis of the influence of strong short-range dis¬ 
order A(po) < 1) was attempted in Ref. [HI via the self- 
consistent Born approximation (SCBA). The latter how¬ 
ever overlooks the diagrams with crossed disorder lines. 
In WSM, these diagrams happen to be of the same or¬ 
der as the diagrams included in the SCBA. Therefore, 
the SCBA is an uncontrolled approximation for a WSM. 
A fully consistent analysis was performed in Refs. IHII^ 
(by means of the ^-expansion). Crucially, it was shown 
that the disorder operator becomes relevant if the initial 
coupling exceeds some critical value 

A(po) > A* - 1. (21) 


For A(po) > A*, the disorder coupling grows with de¬ 
creasing momentum, which leads to the metal-insulator 
transition. If the initial disorder strength is small, 
A(po) A*, then the disorder operator is irrelevant. The 
corresponding running value of the dimensionless cou¬ 
pling constant is 


X{p)=X{po)-, (22) 

Po 

where p ^ maxj^F/'^, T/i;}. We will focus on the case 
A(po) ^ 1 and discard the renormalization of the Fermi 
velocity and field operators. 


2. Long-range disorder 


As one can deduce from Eq. dUD, the long-range 
(po s/v) disorder does not present a problem since the 
corresponding Ej part is much smaller than the 1/{st) 
parameter. Therefore, here we also discard the renormal¬ 
ization of the Fermi velocity and field operators. 


B. The general Dyson equation and its solution in 
the isotropic case 


Surprisingly, even for an anisotropic system, the 
Green’s function preserves its simple isotropic form (in 
the limiting cases of short-range and long-range disorder, 
see Appendix F with detailed calculations). The Green’s 
function has the form^^ 


G^’^is,p) = 


e + f II per 




(23) 
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We will need the following scattering times 


The conductivity then takes the typical Fermi-liquid form 


1 

^impV 

fdn 

T(nno) 


4.TT ' 

1 

^imp^ 

f dfl 

Ttr(nno) 

27r^^p^ J 

Att ' 


_ 1 1 



~ r Ti ’ 


1 

'^impV 

f dn 

Ttr2(nno) 

^TTvjPo J 

47r' 


1 1 

r T2 ' 


^P£(n-n'),no ^nn') 

(24) 


The derivation of the scattering rates is sketched in Ap¬ 
pendix B. In accordance to what was said about the valid¬ 
ity of the Fermi-liquid approach, we only need the value 
of the vertex function at p = . 

This allows us to perform a partial momentum integra¬ 
tion in the Dyson equation and write down the equation 
for the on-shell (p = Ps) vertex (see Appendices A and 
C). It reads 


2^2 


J^^^(nno) =a+ 

q 


47ri;j|pg^2 

dfl' (1 + n'o-)J-^^^(n'no)(l + 


Att 


r(n'no) rtr(n'no) 


(25) 



— OO 


s^ds 

27r 


ds tanh 



(29) 


Now, we write down the temperature dependence of the 
conductivity in the cases of the long- and short-range 
disorder potentials. 


C. Short-range potential 


In the case po ^ max{T/i;, e^/v}^ the structure factor 
^Pe(n-n') depends only slightly on the scattering angle 
(between the n and n' directions) as the scattering is 
nearly isotropic. This means that g{p^{w — n')/po) ~ 
p(0) = 1. The transport scattering time is 


imp'^O 2 


Ttv^{e) 3r(e) 
and the conductivity reads 


37r'i;^PQ 


1 e^i;^Po 


^imp ^0 


(30) 


(31) 


The conductivity is independent of the chemical potential 
and temperature in the short-range disorder limits. 


D. Long-range potential 


where we denoted = J^^^(nno), meaning 

that the angular dependence of J^^^( 5 ,np£) is clearly 
defined just by the angle between the anisotropy axis and 
momentum direction (see Fig. [ 6 ] illustrating the scatter¬ 
ing angles). The energy argument in the scattering times 
is suppressed. In the isotropic case, equation (|25]) allows 
for a simple solution. Since J is a renormalized current 
operator, we look for a solution in the form of the polar 
vector ansatz 


jARR{s, n) = Ji{s)a + J 2 (^)n + Js{s)n{an). (26) 

The ansatz ([26]) turns Dyson equation ([25]) into an al¬ 
gebraic one. Changing the integration measure in the 
expression m for the conductivity p‘^dp = e^dpjv^ ^ we 
arrive at the following suitable formula for the conduc¬ 
tivity tensor 


<^a/3 — 


OO 


tanhJi(e) + J 2 (£) + J3(£) 


27r 


2T 


r(e) ^ ri(e) 


(27) 


The Ji + J 2 + J 3 term is extracted with the help of 
Eqs. (I2K1) and 


Ji + J 2 + J 3 = \^ 


(28) 


This case corresponds to po ^ max{T/'r’||,5 f/'^||}- 
Now, the scattering is strongly anisotropic. Then, 
P£^(n — n')^ = 2 p£^(l — cos 6 >) ^ Pe^0‘^, and 






Stt 


ple‘^ ’ 


(32) 


where gi = f^g(x)xdx (we assume the convergence of 
the corresponding integral). The conductivity then reads 


8 

<^a/3 — - 


4 2 

SpPo 


3 nimpitggi v2 




f(x) = 


.2^2 


IV 


(33) 


In particular, when the doping level is low, sp T (x ::$> 

1 ), 


we have 


f(x) « —X , 


56 e2p2 

^a/3 ^ /I r 2 2 • 

45 nimp'Uo^l 


(34) 


(35) 


The temperature dependence a ^ for long-range dis¬ 
order of arbitrary form can be qualitatively explained 
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in the following manner. The Born transport scattering 
rate is oc {pos)~‘^ and the density of states satisfies 
v (X e^. For T the characteristic energy of the 

charge carriers becomes e ^ T. Therefore, the Drude 
conductivity takes the form 

cr ocrtr(T)z/(T) ocT"^. (36) 

Formula (ESI surprisingly gives the same dependence 
even for Coulomb disorder. As was pointed out in 
Refs. l2in the screening length for the Coulomb dis¬ 
order potential is temperature-dependent and po ^ T. 
The potential amplitude, however, also depends on tem¬ 
perature 1^0 ^ ^ e^po oc T and the dependence 

survives. 

Next, our attention turns to the anisotropic case. As 
we will see, the longitudinal and transverse conductivi¬ 
ties exhibit different temperature dependence. Also, the 
uniaxial anisotropy introduces additional geometric fac¬ 
tors, which are usually needed by experimentalists. In 
some cases, we managed to obtain exact results. 



Figure 6: (Color online) Scattering event: q is the 
momentum of the incoming particle, no is the axis of 
the uniaxial anisotropy. For simplicity, the coordinate 
system is chosen such that q, no span the xz plane; p is 
the scattered momentum, 0' and p are the polar and 
azimuthal scattering angles, while y is the scattering 
angle with respect to the anisotropy axis. 


IV. ANISOTROPIC CASE 

A. The Dyson equation and conductivity 

For convenience, we now introduce a modified ver¬ 
tex function I = J^^^/(2/r — l/rtr) and denote 

= I(nno) 

T( \ ( ‘^ 1 ^ _ ^imp'^o^^ 

U(nno) rtr(nno); “ ^ ^ 

(37) 


B. Short-range potential 

1. 5-correlated potential 

A weak momentum dependence of the potential leaves 
it essentially isotropic, because it does not depend on the 
momentum. Therefore, the answer for the ^-correlated 
potential (its characteristic momentum is po = oc) is 
obtained immediately 

<y\\ =cr, crj_= , (40) 

where the conductivity a is defined in (|3T]) with v = v\\. 


The solution to the Dyson equation is sought in the 
form of the most general polar vector composed of cr, n 
and no: 


I = her + hrv + /3n(ncr) + /4no(ncr) 

+/5n(noo-) + Zeno + /7no(noo-). 


(38) 


Plugging (|38l) into expression ([15]) for the conductivity 
tensor, we obtain 


a(/3(^) 




27r^ 


^11 




)«a«/3 


(39) 


--5'(nno)(no„n;3 + n^no/j) 


where ^ = A -h /2 + Z 3 + T = /4 + Jo + I7X, and 
X = cos^. 


2. Finite-range potential corrections 


Now it is clear that the anisotropy in the rescaled 
basis enters the conductivity expression only as a 
max{5F,T}/po correction. The first correction in the 
max{5F,T}/po series is easy to compute, expanding the 
Dyson equation. The spectrum also acquires anisotropic 
corrections. They are proportional to the disorder po¬ 
tential amplitude and are irrelevant for sufficiently weak 
potential. We expand the disorder form factor © ac¬ 
cording to 


g{x) ^ 1 -h dg{0)x, x = ^ < 1. (41) 

Po 

Using Eq. Ei) and solving the Dyson equation EH) with 
the ansatz (|38l) . we obtain the corresponding correction 
to the conductivity (the details of the computation are 
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Figure 7: (Color online) A part of the Fermi sphere 
(shown in orange). The small dark sphere represents 
the possible change of momentum due to scattering. 


summarized in Appendix E) 


cr\\ = cr 


1 - dg{0) 


cr_L = 


1 - dg{0) 


1 + 4^-^ 7r^TV3 + 4- 
5 {^'||Po)^ -I ’ 

7 + 8g-^ 7r^TV3 + £|- 
15 (^'||Po)^ - ' 


(42) 


It is worthwhile to note that the anisotropy manifests 
itself in the temperature-dependent conductivity ratio 


0-_L 


= 1 - 


±{1 + C^)dg{0) 


TT^rys + 4 

(v\\Po)^ 


(43) 


two coupled differential equations 

$"(1 - x^) + x$'(-3 - y + - 1)) 

_$(1 + ^2 _ ^2(^2 _ 1)) ^ _ 1 (^2 _ (^2 _ i)^2)3/2^ 

^-"(1 - x^) - a;[(y + 1) - - 1)]^' = -2$', 

(45) 

where, for brevity, we switched to dimensionless functions 
^^(^,^') (^,^') and go = g{x)dx. Natu- 

rally, the problem of boundary conditions immediately 
comes on stage. As will be shown below, it is possible 
to circumvent it in a number of important limiting cases. 
The Dyson equation for T determines the T(x) function 
up to an arbitrary constant. This, however, does not 
cause any difficulty for the conductivity, because T(nno) 
enters the integral (|39l) with n as a multiplier and forms 
an odd function of the polar angle. 

The resultant expressions for (t\\ and read 

1 

c^ll(^) = o’i/(F^ j[^{x)x‘^+ ^{x)x\dx, 

-1 

aj_{T) = ^-p-f(^^^ J ^{x){l-x^)dx, 

-1 

8e2e4p2 

2 2 ’ 

where the function f{x) is defined in (|33l) . 

Next, we explore the geometrical implications of 
anisotropy in the three limiting cases. 


Equation (|^ is one of the central results of the paper. 
We have just found that the conductivity components 
in the parallel and transverse directions exhibit different 
temperature behavior. Although Eq. (|43|) is obtained in 
the limit po max{T/u, we guess that the tem¬ 

perature dependence of the cr\\^‘^/a± ratio should become 
more pronounced when po ^ max{T/u, 


C. Long-range potential 

The long-range potential u{p) corresponds to the case 
Po maxj^F/'^ii, T/u||}. Now the ratio 

(44) 

Ps 

becomes an additional small parameter of the problem. 
The disorder potential does not change much the momen¬ 
tum of an incoming particle Sp ^ po ^ Ps- Therefore, 
hereafter we can consider small scattering angles wher¬ 
ever necessary (see Eig. O. 

The calculations are rather cumbersome. They are pre¬ 
sented in Appendix E. The Dyson equation is turned into 


D. Limiting cases of anisotropy 

1. Easy plane, ^ 1 


The first one is the case of strong anisotropy, when 
^ 1 (u_L ^ '^ii). Then, the Dyson equations (|^ reduce 

to the first-order differential equation 


(i>a;)' = I \/l - a;2 , 

T ^ const + O 

The constant in T(x) is irrelevant as was mentioned ear¬ 
lier; while the solution, which is analytic in the interval 
X G [—1,1], reads 



$(x) = |(y 


1 — x‘^ -\- 


arcsm x 


4 \ X 

The conductivity can be then written as 
37r^ ./T' 


)■ 




(48) 


( 49 ) 
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The conductivity ratio becomes 


0-_L 


3 

i+41n2 


O(r'), ^^oo. 


(50) 


2. Weak anisotropy, ~ 1 

Now we turn to the case of v± ~ v\\. Let ^ = 1 + 
^ 1. We expand the vertex function according to 
$ = (the value 1/4 corresponds to the isotropic 

scattering). The corresponding equation reads 

- x^) - Ax5^' - 2(5$ = -^^(1 - x ^). (51) 

The solution analytic in X G [-1, 1] has the form 



ln((T||), 

ln(crj./f^) 


ln(r/ep) 


^(a^) = \- - 5), 

^(^) = -(<5C)|. 

Then, the conductivity reads 

It is worth noting that to the first order in the deviation 
from anisotropy, the ratio 

^ = 1 + 0{5i'^), 6^ < 1. (54) 

does not exhibit any shift. 


(52) 


(53) 


3. Easy axis, ^ ^ 1 

In this case, we are unable to obtain exact coefficients. 
However, it is enough to say that at ^ = 0 there exist 
finite solutions and T(x), which lead to some finite 
integrals with and x in the [—1,1] range. Then, we 
have a qualitative answer 



where ci and C 2 are constants of the order of unity. Then, 
the conductivity ratio takes the form 

^ = £i + e)(^2), ^<1. (56) 

C2 

The above results are summarized in Fig. [51 


Figure 8: (Color online) Conductivities cr||^^(T, ^) as a 
function of the anisotropy parameter tana = 4, 
tan /3 = 1. 


V. CONCLUSIONS 

The results of this work can be summarized as fol¬ 
lows. We have rigorously studied the conductivity of 
isotropic and anisotropic weakly disordered WSM in two 
important limits: short- and long-range disorder (po ^ or 
^ max{T, ^pj/r’ii). With the help of the diagrammatic 
approach, we have been able to justify the applicability of 
the Fermi-liquid theory to WSM. The disorder potential 
was assumed to have a general form. 

We have found that short- and long-range disorder 
leads to different temperature dependences of the con¬ 
ductivity. In the case of short-range disorder, we have 
discovered that uniaxial anisotropy leads to even different 
temperature dependences of the longitudinal and trans¬ 
verse conductivities. In contrast, the long-range disorder 
yields identical temperature dependences of conductivity 
components, ^ T^, for T ^ 

We have also explored the dependence of the conduc¬ 
tivity tensor on the anisotropy parameter ^ and estab¬ 
lished general scaling relations in the cases of strong and 
weak anisotropy. We have managed to compute analyti¬ 
cally the geometric factors for the conductivity tensor in 
the limit of strong ^ ^ and weak ^ ~ I anisotropy. 

The recent experimental data on WSM makes it pos¬ 
sible to estimate the Fermi velocities and Fermi ener¬ 
gies for typical samples. For example, for both NasBil 
and Cd3As24, the Fermi velocity ratio is ^ ~ 4; thus the 
regime ^ ^ I is realized and the anisotropy of the conduc¬ 
tivity should be clearly pronounced (^ <^^). The Fermi 
energy takes values in the lOO-IOOO K range. Therefore, 
depending on the disorder correlation length, results (H9l) 
or (|l3|) should be applicable. 

The developed approach and the obtained results pro¬ 
vide a good basis for further progress in the field of trans- 
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Figure 9: Contour for polarization operator 11 (cj). 

port phenomena in 3D systems with Dirac points in their 
energy spectrum. 
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Appendix A: Expression for the polarization 
operator 

In order to work out the polarization-operator diagram 
in Eig.[3](left or right), we follow the scheme proposed by 
Eliashberg^. To build an analytical continuation of the 
expression for the Matsubara polarization operator, we 
need to establish the analytical properties of the vertex 
function Ja(^, z-\-iuJn^ iojn] p) in the whole domain of the 
complex variable This is done via the Lehman repre¬ 
sentation and it was discussed in detail in, e.g.. Ref. [^. 

The conclusion is that the domain of analyticity of 
Ja{z,z -h iujn^iujn^p) is a complex plane with two hor¬ 
izontal cuts: Im( 2 : -h iuj) = 0 and Im(z) = 0. Since we 
need a retarded vertex function, we put > 0. Next, 
the three vertex functions are defined in accordance with 
the structure of the cuts 

z -\-if Im2: > 0, 
z if — icjn < Im2: < 0, (Al) 
z-\-iuj^iuj) if Imz <—iujn- 

The general expression for then becomes 


na/3(iw„) = -T '^Ja{iek,iek + p)G(iejk + iojn,p)jpG{iek,p) 

f ds E 


(A2) 


The contour C is shown in Fig. 0 As usual, the integral over the large circle vanishes and we are left with integrals 
over different branches 

na/3(iw„) = J ^tanh^|jf^^(£,£ + iw„,iw„;p)G'^(e + iw„,p)j>G^(e,p)- 

- J^^^{e,e + iujn,icon;p)G^{e + iLo„,p)j/3G^{e,p) + - iLo,e,iujn,p)G^{e,p)jpG^{e - iw^.p)- (^3) 

- iujn, s, ioJn, p)G'^(e, p)jpG^{£ - i<^n, p)}- 
Making the analytic continuation iujn u iO, we obtain 


nf/jH =tr J ^^^i{G^(e,P)-^a'^'^(e,e + w,w;p)G^(£ + w,p)cT/3 


, E + UJ ^ E 

tanh ——-tanh —- 

2T 2TJ 


+G^(5, p)J^^^{e, 5 + cj, cj; p)G^{e + uj; p)cr /3 tanh — - G^(5, p) e-\-uj,uj; p)G^(5 + cj, p)cr /3 tanh 


(A4) 

E UJ' 


'■y 


Here, the vertex functions are defined diagrammatically in Eig. [3t. They obey the Dyson equation presented in 
a diagrammatic form in Eig. [H The disorder is static and does not cause a change of frequency in the diagrammatic 
loops. Thus, the Dyson equation takes an especially elegant form 


T Q) — T 


Po 


ij 


dp , 

S^q-p.no 


(27r 


G(isk, p)Ja{isk, isk + iun, q)G{isk + p) • (A5) 
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We are interested in the zero-frequency response of the system and set the external frequency cj = 0. After the 
analytic continuation, the vertex function is split into RRR, AAR and ARR parts 


=(Ta + 
= (T„ + 
J^^^(e,q) = CTa + 




j P)Ja^^{£, P)G^{e, P) , 




/ j^9'q-p,noG^{e,p)J^^^{e,p)G^{e,p), 




dp , 




(A6) 


The learned reader already knows that only 
(so called singular) vertex undergoes a strong (non- 
perturbative) renormalization due to the interaction with 
impurities, while the other two exhibit weak perturbative 
corrections. The reason is, of course, the position of the 
poles of the Green’s functions. Since we claim some rigor, 
we discuss this issue in detail in Appendix B. 


Appendix B: One-loop structure of the Green’s 
function 


1. General expression 


The general expression for the one-loop self-energy 
reads 

yR, . ^ riimpul f dp s + v\\pa 

^ J (27r)3 (e + zO)2-i;2p2yq-P,no 

(Bl) 

and its imaginary part 

■2-2 r ^ 


Iml]^(e,q) = — 






(B2) 


The self energy has the following tensor form 


q) = a(£, q) + /3(e, q)cr - 


icrs(5, q) 


2r(e,q) 2Ti{e,q) ’ 
(B3) 

where s and (3 are vectors in the (q, no) plane and 

am /Pon||N2 min{pon||,£} 1 ^ ^ 

s' e \ e ) e er 

For the case of isotropic potential, the anisotropic short 
range potential (po ^ Pe) and anisotropic long range 
potential (po ^ Ps)^ vector s = n and the scattering 
rates r and ri are immediately extracted and presented 
in (EH) . 

As is seen from the structure of the integrals m and 


T{e,-q,0) = T{e,q,0); 
s{s,-q,0) ^ s{s,q,0) 

Ti{e,-q,0) Ti{e,q,0)' 


(B5) 


(a) 


(c) 


G^{s,q} ■ 


G\e,q)G\e,q) ' 




(b) 

£ 

—s 

s 

—£ 



G^(e,q)G^(e,q) ' 





Figure 10: Poles of (a) G^{e,q), (b) G^{e,q)G^{e,q), 
and (c) G^{e,q)G^{e,q), in the q plane. 


Here 0 is the angle between q and no and the parity 
relation (|B5D follows from the measure invariance under 
the change 0 ^ tt — 0 ^ (p ^ (p tt. 

Now we may write the Green’s function in the standard 
way 


G«{£,q) 


e - a + ^ + cr(q + f3 - 
(e-a + ^)-(q + ^-^) 2 - 


(B6) 


In a typical computation, we need to integrate prod¬ 
ucts of the type G^{s,q)G^{s,q), G^{e, q)G^{e, q) and 
Gme,q)G^{e,q) over the momentum q G (—oo,oc). 
Therefore, we need to know the position of the poles of 
the Green’s function in the g^-domain. 

The retarded function G^(5,q) obeys the causality 
condition, which means that its poles, defined by the 
equation 


det[5 — v\\qcr — q)] = 0 , (B7) 

lie in the lower half-plane of the complex variable 5, i.e. 


rvsie,e,0) 

+ ; A > 0 


r{e,e,0) Ti{e,e,0) 


(B8) 


for any 5 and 0; n stands for q_/q. Since we need the 
position of the poles in the g^-plane, we solve the spectral 
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equation (|B7D for q and obtain two roots. Their imagi¬ 
nary parts are 


lmq{€) = ± 


1 

-2r{s, ±s, 0). 


s(5, ±5, 0)n 

2 ri(e,±e,0) ‘ 


(B9) 


Using ()B5D and ()B8D . we see that the position of the 
poles of G^{s,q) in the complex plane q is always such 
as depicted in Fig. [TOh . The poles of G^{e,q)G^{s,q) 
and G^{e,q)G^{e,q) are sketched in Figs. [TOb andfTOb. 
Now, we see that while integrating the product G^G^ or 
G^G^ over q the contour of integration can be deformed 
to pass far from the poles. When integrating the product 
the contour is squeezed between the nearby poles 
and the deformation is no longer possible, leading to their 
enhanced contribution. 


at p = s/v\\ as well as explore its asymptotic behavior 
at p ^ oc (to justify the Fermi-liquid approach). The 
asymptotic behavior is easier to extract focusing now 
on the isotropic potential. This allows us to capture all 
the principal details and avoid unnecessary complications 
due to uniaxial anisotropy. 

The ansatz for the vertex comprises only three vertices 
(Eq. 1^61) . However, in formula (f^6l) the vertex is taken 
on a mass shell at p = p^ = ^/u. Here we consider a 
complete vertex 


J(e, q) = Ji{£, q)a + J^ie, q)^ + Jz{e, q)^^2 


(Cl) 


Appendix C: Dyson equation 


1. Structure of J"^^^(£,p) 

In order to find the conductivity, we need to solve 
the Dyson equation for the vertex function J^^^(5,p) 


For brevity, we will omit the energy symbol 5 in J{e^ q) = 
J{q) to restore it later. Using the identity [e-\-VY>cr](T[e^ 
upcr] = {e^ — v‘^p^)(T + 2up[5 + uptr], we obtain the equa¬ 
tion 




q(qo-) 


J^{q)cT + Uq)^ + JM „2 

dp 


= cr 


pI 


j 


q" q- 

Mp) + 2p(Pe + per)] + J 2 (p)^ + 2p£(pcr)] + Jsip)^ {pe'^+p‘^){pa)+2psP^ 


(27r)' 


rdq-p 


(yp ri(e,p)) _ r(£,p)) n (e,p)) 


(C2) 


For the integral (|C2p to converge, the vertices Ji, 2 , 3 (p) 
must be bounded when p ^ oc. As in expression for the 
conductivity, the integrand has a sharp peak at p = p^ 
of width Ap ^1/ (vt) near the Fermi surface. Generally, 
we cannot restrict the computation of the integral by 
expanding the integral near the Fermi surface. Such an 
expansion leads to a significantly deformed integrand (see 
Fig. [TTT aU. First, however we discuss the vertex at p = 
Ps- 


2. On-shell vertex: 

If we are interested in the value of Ji, 2,3 at p = p^ 
then the range of momentum integration is the sphere of 
convergence of the disorder structure factor p, namely 
the sphere of radius po with the center at the exter¬ 
nal momentum q = p^n. There are two distinctive 
contributions. The first one comes from the integra¬ 
tion inside the small sphere: |p — np^j < l/(ru), where 


the integrand is peaked due to the proximity of the 
poles of G^G^; the second one comes from a thick shell 
|l/(TUf)| < |p—p^nl <po. Alternatively, one can think of 
the first contribution as the one coming from the singu¬ 
larities of G^G^ and the second of from the singularities 
of gnpe-p' Simple estimates give 

J dpG^{p)JaG^{p)gnp,-p ~ Te^Ja{Pe)gp,, 

\p-npe\<:k 

J dpG^{p)JaG^{p)gnp^-p ~ 

(C3) 

Comparing the two contributions, we see that if the dis¬ 
order cut-off Po is not too high, the contribution from the 
thick shell is I/{re) x (po/Pe)^ smaller than the one from 
the small sphere due to the G^G"^product. Therefore, 
one can indeed discard the contribution from the shell 
^ \p ~ Ps'^l ^ \Po\ (or equivalently, the contri- 
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Figure 11 : (Color online) (a) Approximate integrand 
G^{p^^)G^{p)Ps‘^9pein-n') (transparent) versus the 
exact one G^(p, 6 >)G"^(p, 6 >) 5 fq_pP^(solid), for a typical 
potential, (b) The range of integration over p. The 
transparent region (radius po) corresponds to the 
smaller contribution. 


but ion from the poles of the potential gnp^-p)' Thus, in 
the leading order of the l/(T5)-expansion, one can take 
into account only the poles of the G^G^ term. 


3. Vertex at p ^ oo 

The same logic of estimates holds for the large p- 
behavior of the vertex. We now plug ansatz (PTj) into 
the Dyson equation ()A 6 p and employ the estimates (|C3D . 
modified for the external momentum q ^ Pe, po- The 
computation is simple though tedious. Solving the sys¬ 
tem of self-consistent equations, one arrives at the fol¬ 
lowing asymptotics 

Ji{q) = l + l3igg^J, 

Po 

2 

Mq) = P2—Podqgq^J, (C4) 

Ps Po 

Mq) = ^3gq^J, 

Po 

where J = Ji{ps) + J 2 {Pe) + J^{Pe) and / 3 i, 2,3 are con¬ 
stants of the order of unity. We see that the estimate for 


the vertex in the isotropic case can be rewritten as 

J(q) - (T ~ rfi((g), (C5) 

where r is some restricted vector function of the momen¬ 
tum q. Therefore, the renormalized current vertex J(q) 
approaches the bare one cr at large momenta q with the 
rate proportional to the structure factor g{q). 


4. On shell Dyson equation 

Now we are able to write down the Dyson equation in 
a somewhat simplified form performing integration over 
the absolute value of the momentum. As was argued pre¬ 
viously, the most important contribution comes from the 
peak in the G^G^ product at p' = p^. Therefore, we 
put p' in the argument of J( 5 ,p'n') and p'(p — p'n',po) 
inside the integral equation (|A 6 D equal to p^. More¬ 
over, we take into account the fact that in all the lim¬ 
iting cases considered (isotropy, anisotropy with short 
and long-range disorder) the Green’s function assumes a 
simplified form (P3|) rather than (|B 6 p (see Appendix F). 
Changing p^ = C, dp = v\\d(/2e, and changing the lower 
limit of integration to — oo, we perform the integration 

CXD 


p^dp 


■ / . \ 2 / • \ 


^2 ^ \ 2n 



) ] 

TT 

1 

(C 6 ) 

II r 



This way we arrive at the Dyson equation (ESI) and dSil) 
for the modified vertex. 


Appendix D: Estimate of the regular term in 
conductivity 


Thanks to the factor — p^ in the nominator of m 
we can substitute 


(p-£^)G^(e,p) fa (e + po-) l + e)(^)^ 


(Dl) 


and the contour of integration is no longer squeezed be¬ 
tween the proximate poles of the product G^G^. The 
path therefore can be deformed in any suitable way. Es¬ 
timating at p ^ max{p£,po}, G^G^ ^ 1 /^^, we obtain 


I 6ag\max{ps\j{s,ps) - o-|,po|J(5,po) 


<r|}l (D2) 


On the other hand, the convergent expression II is com¬ 
pletely determined by the behavior of J^^^( 5 ,p) in the 
vicinity of the peak and one can substitute Ja{s,p) with 
Ta(^F,PF)- The integral is then estimated as 


II ~ (5„;3|j(e,Pe) - (t|^. (D3) 

Combining Eqs. (|D3D and (IC5p . we obtain relation (HZI. 
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Appendix E: Conductivity, short-range potential ing averages 


In the (^-correlated (momentum-independent) case, the 
vertex functions 4^ and are angle-independent and, as 
follows from Eq. ([^ 


4>(nno) = <ho = Y’ ^(nno) = 0. (El) 

In the case of small momentum-dependent correction 
to form factor dm), these also acquire small angle- 
dependent corrections 

<l>(nno) = 4>o + (54>(nno), 4^(nno) = ^4>(nno) (E2) 


^ _dgifi)p^ 


f dfl' , 

J 47^ 5'p£(n-n'),no'^ 


[(1 -C ^)(nno)no +C ^n] , 


^imp^ 

^ ( ^al3 


fdfl', 1 

(l-r')(nno)2 + 


1 + 


dg{0)pe^( 2 n / ^2 , 12 + 9 ^ ^ 




+ ^)("“o)noano/3y 


(FI) 


of the order of (pe/po)^- Thus, the formula for the con¬ 
ductivity (see Eq. ®) becomes 


Therefore, we see that the self-energy retains its isotropic 
form with the corrections to the scattering times: 


^a/3 — ^ 

/ (54>o — (^4>2 J. 3(54>2 — (^4>o + 25^1 \ 

^ H ^a/3 H ^0q;^0/3 J 5 

(E3) 

where = f^-^{S'^^6^)x'^dx. Now (54>^ and 

(5T^ are easily extracted perturbatively from the Dyson 
equation (see the next Appendix for details). The cor¬ 
rections can be written as 


6a±_ = 


6a\\ = 


2 — ^4>2 

^^2 + 


2r 


(E4) 


T T Pq \ L 3 


2e 


-2 


(F2) 


,1 l%(0)p,V5,, ^-2^. ^2 , l + 7r' 


Eollowing the same steps as in the isotropic case, we ar¬ 
rive at the following equations for the vertex functions 


5^x) = -‘^dg{0)^{l + 4r") + 


Po 

s^{x) = -^%(o)^(i - 

10 Pq 


S'^o + (5^*1 

X H- X . 


(F3) 


We thus obtain the conductivity as a function of the en¬ 
ergy 


5(j±_ = 
(5(711 = 


5 Po g 


(E5) 


Plugging these into da, we determine the temperature 
dependence of the conductivity (|42]) . 


Solving Eqs. (|E3D in a self-consistent way, we find 

= -^r(l + AC‘^)dg{Q)^, 

5 Po 

<5$2 = -^r(l+4r')a5(0)^, (E4) 

5 Po 

<5^'i = 

iO Po 

Erom the last equations, we easily recover Eqs. (|E5E 


Appendix F: Solution of the anisotropic Dyson 
equation 


2. Long-range disorder 

Expanding in the scattering angle, we have 


1. Short-range potential 


We now need to compute the anisotropic corrections 
to the spherically-symmetric self-energy. Erom the struc¬ 
ture of ()B2E we see that we need to compute the follow- 


(n — n')no = (1 — cos O') cos 0 -f sin 0' sin 0 cos p 
^ 0' sin ^ cos p, 

(n — n')^ = 2 — 2cos6> ^ 0'^, 

do' = siiiO'dO'dLp ^ O'dO'dif . 


(E5) 
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Next, we extend the upper limit of integration over 0' to 
+ 00 , and we obtain 


+ (JxO' cos ip\ 


ImS(£,pen) = 6^2 [ ^ l'^ + <^y^'sin<p 

4nvlpteJ 27r ^ 

f°^0'de' ,_2, .2. 2 .- 2 A 

X 7 


crn 


2T(nno) 


( ----- ) • 

VT(nno) Ttr(nno)/ 


(F6) 


In addition to (l2^ , we will need one more scattering time 


1 

Ttroie) 


^imp'^ '^0 
2<2^jjpg^2 


f dVt , n \ 2 

J ^%(n-n').no(l-COS0nn')cOS ifi 


(F7) 


(see Fig. [6] for the definition of the angle ip). To obtain 
the explicit formulae for all scattering times, we expand 
the argument of no integration domain dVt. 

Introducing the notation 

1 _ 2 g^^impg0M0 

To{e) ^ Swvf^pl ’ (F8) 

A{e) = A(nno) = [{p - 1) sin^ 9 + , 


to the anisotropy axis. Now we plug the Green’s func¬ 
tions ([23|) and ansatz (|38l) into the Dyson equation (|37l) . 
The most obvious step in solving the Dyson equation (|37j) 
is to take the zero order in the /^-expansion and to set 
n' = n in the peaked integrand. One immediately ar¬ 
rives at an algebraic equation for all the constants. De¬ 
spite the fact that one is able to immediately obtain some 
of the vertices: /i = 1, /s = /y = 0, we also obtain a 
meaningless equation I 2 = I 3 = 1 12 and the vertices 

/ 4 , Iq remain undefined. This is a signature of the fact 
that one needs to take the next term in the /^-expansion. 
The K expansion in turn leads to the emergence of the 
derivatives of functions I in the structure of the equation. 

We use the commutation relation and find (1 + 
na)a{l + ntr) = 2n(l + ntr). 

(1 + n'o-)I(l + n'a) 

= 2(/i + /2 + /s + / 5 n'no)n(l + n'a) (FIO) 

+ 2(/4 -|— \-Iq + /7n^no)no(l + n^cr). 

Whenever necessary, we expand 

n' = n + (in, 

/(n'no) = /(nno + ^nno) 

= J(nno) + no(5n5„„o/(nno) + dl^J{nno) + ... 

(Fll) 


where go is defined after Eq. ( 03 , we obtain results for 

the cross-sections Next, we perform the integration over angles 
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(Here gi is defined after Eq. (1321) 1. As one can see from 
the expression for Im E (|F6l) , the Green’s function retains 
its simple structure (|^ . One needs to pay particular at¬ 
tention to the fact that the scattering time is now a func¬ 
tion of the direction of the incoming particle with respect 
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The Dyson equation is split into seven differential equa¬ 
tions. In terms of ^{x) and T(x), we rewrite it as follows 
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Summing up the equations, we obtain 
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